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The electronic structure and electron-phonon coupling in Y2C3 is investigated using density 
functional calculations. We find that the Fermi level falls in a manifold of mixed character derived 
from Y d states and antibonding states associated with the C dimers in the structure. Calculations 
of the electron-phonon coupling for Y and C modes show that the former provide most of the 
coupling. Modes associated with C-C bond stretching have large matrix elements, but make small 
contributions to the coupling because of their high phonon frequencies. Substantial electron doping 
of the C-C antibonding states would yield a large increase in the coupling and critical temperature, 
perhaps to values comparable to MgB2. However, it seems unlikely that a modification of Y2C3 
with much higher filling of the C-C antibonding states can be stabilized. 



The rare earth sesquicarbides, ^Cs, where R is a rare 
earth or Y are a family of superconducting materials that 
form in the bcc PU2C3 structure under high tempera- 
ture conditions. Samples produced by arc melting have 
critical temperatures, T c , of 10K to 12K, with little de- 
pendence on the particular rare earth. For example, as 
synthesized samples with the largest and smallest rare 
earths, specifically, La2C3 and Y2C3 have critical tem- 
peratures of 11K and 10.5K, respectively. 1-5 However, it 
was found early on that heavy Th doping raises T c , up 
to 17K in the case of the Y compound. 6 Very recently, 
Amano and co-workers 7 reported that T c for Y2C3 could 
be increased to 18K by synthesis under high pressure (4- 
5.5 GPa). The resulting samples have the same PU2C3 
structure as previously obtained samples but a slightly 
different lattice parameter, suggesting that the effect of 
pressure is to change the composition of the samples. 
Furthermore, they reported that the results were sensi- 
tive to the sintering conditions. The lattice parameters 
varied in the range a— 8.181A - 8.226A, with reported T c 
spanning the range 15K to 18K. The extrapolated upper 
critical field of the 18K samples was estimated to be ~ 
30T. 8 

The PU2C3 structure has four formula units per prim- 
itive unit cell, spacegroup /43<i, Y on site 16c (u,u,u) 
and C on site 24d (v, 0,1/4), where the coordinates are 
in terms of the bcc lattice vectors (see Fig. 1). As may 
be seen, C dimers with short bond lengths are a charac- 
teristic of this structure. Electronic structure studies of 
other metals with C dimers both at the tight binding 9 
and density functional level 10 suggest that the electronic 
structure could have substantial C-C antibonding char- 
acter near the Fermi level. Such states may be expected 
to have extremely high deformation potentials, reflecting 
these very strong C-C triple bonds. Also reflecting the 
very strong covalent bonds, phonons that modulate the 
C-C distances would be very stiff, which, if these states 
were responsible for the superconductivity, would yield a 
moderate electron-phonon coupling, A, and a very high 
logarithmic average phonon frequency and therefore high 



prefactor in the McMillan equation for T c . In this sce- 
nario, a very strong dependence of T c on stoichiometry 
(more precisely, occupation of the antibonding states) 
might be anticipated, as well as the possibility of very 
high T c in more optimally doped samples, perhaps with 
values like in MgE>2 and ^Ceo- Here we report density 
functional calculations of the electronic structure and the 
coupling constants of a C-C and an Y mode aimed at de- 
termining the extent to which this scenario applies in 
Y 2 C 3 . 

The present calculations were done using the gen- 
eral potential linearized augmented planewave (LAPW) 
method, 12 with local orbital extensions 13 to relax lin- 
earization errors and to treat the semi-core states of Y. 
LAPW sphere radii of 2.0 ao and 1.24 ao were used for 
Y and C, respectively. Well converged basis sets (conver- 
gence tests were done) consisting of 4250 basis functions 
for the 20 atom primitive cell were used, along with con- 
verged Brillouin zone samplings. The core states were 
treated relativistically and valence states semirelativis- 
tically. For the crystal structure we used the reported 
bcc 743d PU2C3 structure with the experimentally de- 
termined lattice parameter, a = 8.226 A, from Ref. 7 
but relaxed the internal coordinates using the calculated 
atomic forces. We found significant changes in the C-C 
dimer bond length from the available experimental refine- 
ment (Ref. 11). Specifically, we obtain a C-C distance of 
1.33 A, and a Y-C distance of 2.51 A, corresponding to 
internal parameters, w=0.0504 and ^=0.2940. The cor- 
responding full symmetry Raman modes are an Y domi- 
nated mode at 175 cm -1 and an almost pure C-C bond 
stretching mode at 1442 cm -1 . It seems likely that the 
difference between our results and the C position of Ref. 
11 is related to difficulty in the refinement due to sample 
stoichiometry, which is an issue in these sesquicarbides, 
as mentioned above. 

The calculated electronic band structure and corre- 
sponding electronic density of states (DOS) are given in 
Figs. 2 and 3, respectively. The Fermi surfaces are dis- 
played in Fig. 5. The band structure agrees qualitatively 
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with a very recent report by Shein and Ivanovskii, 14 al- 
though there are noticable differences in detail, presum- 
ably related to the different C-C bond length in that 
calculation. 

The band structure shows a manifold of 12 C 2s de- 
rived bands extending from -14.9 eV to -6.2 eV, relative 
to the Fermi energy, Ep (note that there are 12 C atoms 
in the primitive cell). This is followed by a narrower 
manifold of 18 C 2p bonding bands, associated with the 
C dimers. This, in turn, is separated by a ~ 2.5 eV gap 
from a broad manifold of mixed C-C antibonding and Y 
character. The Fermi level lies ~ 1 eV from the bottom of 
this manifold. Fig. 4 shows the contribution from states 
near Ep to the charge density. The plot is for the energy 
range that contributes 1 e per unit cell going down from 
Ep- This electronic structure is consistent with a con- 
ventional picture 9 of strongly covalently bonded C dimers 
embedded in a solid that has metallic cohesion involving 
mixture of Y and antibonding C-C states. 

As may be seen, there are several bands crossing the 
Ep. This leads to a complex multisheet Fermi surface 
(Fig. 5), with no obvious strong nesting features. 15 The 
density of states at E F is N(E F ) = 1.88 eV -1 on a per 
formula unit both spins basis, corresponding to a bare 
linear specific heat coefficient, 7& a re = 4.4 mJ/mol K 2 . 
This is approximately 60% larger than the experimental 
value of 7 = 2.8 mJ/mol K 2 (ref. 2), but this comparison 
should be interpreted with caution since the experimen- 
tal stoichiometry was not well established. We note that 
Th doped Y 2 C 3 with T C =17K (closer to the 18K reported 
for Y2C3 by Amano and co-workers), has a reported ex- 
perimental 7 = 4.7 mJ/mol K 2 , and that the DOS is 
peaked at the stoichiometric composition, implying that 
off-stoichiometric samples would likely have lower 7, both 
because of lower jbare and presumably lower electron 
phonon enhancement. 

Considering the apparent absence of strong nesting, 
and the large unit cell (four formula units, 20 atoms), 
it is expected that the zone center modes may be at 
least semi-quantitatively representative for determining 
the electron-phonon coupling. As a further simplifica- 
tion, we focus on two modes in order to obtain the gross 
features. These are the two full symmetry modes related 
to the internal structural parameters. As mentioned, 
these are an Y dominated mode at 175 cm -1 and an 
almost pure C-C bond stretching mode at 1442 cm -1 . 
We determined the mode A for these two modes using a 
frozen phonon method, in which LDA calculations were 
done as a function of the phonon distortion, and the de- 
formation potentials calculated directly from the change 
in electronic structure on the Fermi surface. The result 
is shown as a function of rigid band Fermi energy in Fig. 
6. For this plot, it was assumed that the lower frequency 
Y mode (which modulates both Y distances and Y-C dis- 
tances) is representative of all modes except the 6 C-C 
bond stretching modes, which were assumed to be repre- 
sented by the high frequency C-C mode. Most likely, this 
is reasonably good for the bond stretching modes but un- 



derestimates the phonon frequency for the other modes, 
since the C dimers are lighter than an Y atom. However, 
considering the crudeness of representing the system with 
two full symmetry zone center modes, it hardly seems 
justified to attempt corrections, and instead we view the 
result as semi-quantitative. 

The calculated electron-phonon coupling for the stoi- 
chiometric band filling is A = 0.6 and is dominated by 
phonons other than the C-C bond stretching; the C-C 
bond stretching provides less than 10% of the total A. In- 
serting A = 0.6 into the McMillan formula, with /i*=0.1 
and cj=175 cm -1 , we obtain T C =5K, which is fair agree- 
ment considering the crudeness of the approximation em- 
ployed. 

What is significant is that A has a maximum at the sto- 
ichiometric band filling. This is related to the fact that 
the Fermi energy falls on a peak in the DOS. The depen- 
dence of A on band filling and the difficult synthesis of 
stoichiometric sesquicarbides provides an explanation for 
the experimentally observed variability of T c depending 
on preparation conditions. Going to lower band filling, A 
decreases and then increases again to a value exceeding 
A = 1, reflecting the structure of the DOS, but this only 
happens quite far down, corresponding to an electron de- 
ficiency of about 1 e per formula unit. The C-C strecthing 
contribution is small throughout this region. Going to 
higher electron count, A again decreases initially and then 
increases, again to values exceeding A = 1. However, in 
this case the C-C bond stretching contribution becomes 
significant and as a result the logarithmic average phonon 
frequency also sharply increases. In principle, this would 
mean that heavily electron doped Y2C3 would become a 
high temperature superconductor, like MgB2 and AsCqo- 
However, this is based on a rigid band filling, which is a 
highly unlikely scenario. Instead, population of C-C anti- 
bonding states should be highly disfavored, and therefore 
it would be rather difficult to chemically reach the region 
where this contribution to A becomes large. 

In summary, we report electronic structure and mode 
electron-phonon calculations for Y2C3. The results are 
consistent with conventional electron-phonon mediated 
superconductivity related to Y - C phonons. The results 
show that the superconductivity is qualitatively like in 
other metal carbide superconductors, and not like MgE>2 
or A 3 C 6 o. 

Work at the Naval Research Laboratory is supported 
by the Office of the Naval Research. The DoD-AE code 
was used for some calculations. 
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[52 atoms] [l8-bonds]Y2C3 

FIG. 1. Crystal structure of Y2C3 showing the C2 dimers. 
Large spheres denote Y and small spheres denote C. 
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FIG. 2. Calculated band structure of Y2C3 with the re- 
laxed atomic positions. The lowest bands shown are part of 
the C 2s manifold. The bottom of this manifold is at -14.9 
eV. The lower panel is a blow-up around Ef, which is at 
eV. 
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FIG. 3. Electronic density of states of Y2C3 with the re- 
laxed atomic positions. The solid line is the total density of 
states on a per formula unit basis and the dashed line is the 
Y d contribution as defined by projection onto the Y LAPW 
sphere of radius 2.0 clq. 
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FIG. 4. This figure shows the charge density originating 
from the valence bands near Ef (see text) in 1/4 of the primi- 
tive unit cell around C (dumbells, red online) and Y (intercon- 
nected network, blue online). Isodensity surfaces correspond 
to p(r) =0.00112 e/al for C and 0.00092 e/ag for Y. Small 
differences between the charge density values were needed to 
improve visibility of the two surfaces. For the same puprose, 
the parts of the Y charge density surfaces, penetrating from 
the neighboring unit cells, were removed. Note the metallic 
like Y derived density and the C antibonding density. 





FIG. 5. LDA Fermi surfaces of Y 2 C 3 . 
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FIG. 6. Mode A and average phonon frequency (see text) 
as a function of energy. The heavy dashed (heavy solid) line 
denotes the full symmetry Y (C) derived mode A as a func- 
tion of energy, normalized as if these modes are representa- 
tive. The light dashed line is the total A with this assumption. 
The light dotted line is the logarithmic average of these two 
frequencies with A, which would enter the prefactor of the 
McMillan equation. 
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